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Abstract 

We obtain explicit expressions for the rj and rf masses and decay constants using Ul(3) ® Ur(3) 
chiral perturbation theory at next-to-leading order in a combined expansion in p 2 and 1/N C . A 
numerical fit of the parameters appearing to this order of the expansion is also discussed. 
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The masses of the 1=0 pseudoscalar 77, rf particles as well as the value of their mixing angle have 
long been the subject of discussion and theoretical interest, also in recent years. Phenomenologi- 
cally, the situation of the rj — rj' mixing remains not completely settled, the reason beeing the high 
sensitivity of the mixing angle 9 on the deviation (A) of the rj^ mass from the Gell Mann-Okubo 
(A = 0) relation gj 

m 2 m = ±(4m 2 K -ml)(l + A). (1) 

Indeed, if A = the singlet-octet mixing yields 9 ~ —10°, whereas data on J/ty — * rj{rf) r ) and 
77(77', 7T ) — ► 77 favour a higher mixing 9 ~ —20° a factor of two bigger, that can be accomodated 
from ([l]) with just a small A ~ 0.16, (see also ref. ||, which advocates for a smaller mixing). 

The value of the mixing angle is linked to the value of the 77' mass, which is heavier than the 
octet of light pseudoscalar s, = 958 MeV ~ 2Mk, but still lighter than the next 1=0 candidate 
r/(1279), which is so distant from the 77 mass that it might not significantly mix with it. Therefore, 
we only consider 77 — 77' mixing. 

Within the SU(3)l <8> SU(3)r Chiral Lagrangian H analysis, the shift downwards of the 77 mass 
due to this mixing is an 0(p A ) effect of SU(3)r+l breaking, proportional to (m s — mu + m <i ) 2 anc j 

multiplied by the constant L^f U ^ , which contains information of the 77' state that has been integrated 
out. 

In this article we re-analyse these issues on the basis of large- N c Chiral Perturbation Theory 
(xPT). The 77' is regarded as the ninth goldstone boson of a nonet, together with the octet of 
light pseudoscalar s, and the corrections are treated in a double expansion, both in powers of quark 
masses m quar k and 1/N C . Such a formulation for this problem is plausible because M 2 , — M c 2 ctei is 
an effect induced by the Ua{^-) anomaly, which vanishes in the large- N c limit as 1/N C . Furthermore, 
the fact that phenomenollogically it is found that Mw — M oc t e t ~ M oc t e t, together with the usual 
bookkeeping of quark masses as m quar k ~ 0(M 2 ctet ), suggests that the relative magnitude of the 
double expansion may be regarded as 

m q ~l/N c ~M* ctet ~p 2 ~ (S); (2) 

this is the approach that we adopt and henceforth we shall refer to it as the combined expansion. 
It has already been used in the literature ||], in particular Leutwyler has forcefully pursued the 
analysis of the masses of 77 and 77' in order to discard the possibility that m u = 0. In this note we 
reproduce Leutwyler's results and proceed further to obtain a fit for the masses and decay constants 
of the goldstone boson nonet Q. 

Let us start from the Ul(3) (g) Ur(3) chiral lagrangian which has been recently proposed to 
0(p 4 ) 0, from which we share the notation and conventions, and which we refer to for references 
therein. The nonet of fields are gathered in a unitary 3x3 matrix U S U(3), whose determinant 
det U = exp(7\/6 ? ? //) differs from unity by the presence of the 77 field. The chiral lagrangian 
£( p 0) + £( p 2) + ... reads 

V) = -Wo(x) , 

£ (p2) = W^XXD^D^U) + W 2 (X)(U*x + X ] U) + iW 3 (X)(U^ X ~ X ] U) 
+ W 5 {X){U\D^U)) {ia^). 

r)o, r]g are the eigenstates of SU(3)r+l- Although in the text we are occasionally loose when referring to rf , r\ 
instead of 7/0, 7/s, the distinction is carefully made in the calculation. 
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(3) 



We have only kept the axial source because we take the divergence of the axial current as the 
interpolating fields for the nonet. The axial source appears in the last term as well as through 
the covariant derivatives of the fields: 

The quark masses appear in 

X = 2B diag(m u , m d , m s ) . (4) 

As a consequence of the Ua(1) anomaly, the lagrangian has coefficient functions that depend on 
X = log(det U) = i\^6r] /f, which are not fixed by symmetry arguments. From the large- iV c 
perspective however, only the lowest powers of X survive, for each new power of X is suppressed 
by a factor of 1/N C . Of interest to us, the only terms that prevail are: 

f 2 f 2 
W {X) = Constant + J —v Q2 X 2 + ... , W X {X) = J — + ... , 

W 2 (X) = L + , W 3 (X) = -i^v 31 X + ... , 



W 5 (X) = £ V50 



(5) 



with vq2, V31, V50 ~ 0(1/N C ). We shall expand the chiral lagrangian and keep terms up to 0(5 2 ), 
according to the combined power counting (|). Recall that f 2 ~ 0(N C ), B ~ O(l). Within this 
approximation the chiral logarithms are suppressed (they appear at 0(<5 3 )), although some terms 
of the same order 0(p 4 ) have to be included at tree level. More precisely, the contributions 

£ {p4) = L 5 {D^ D»U (rf X + X ] U)) + LttfUrfU + rfx^X) (6) 

should be included, since L 5 and Lg|| are 0(N C ). Notice that a term L 7 (U^x ~ X ] U) 2 appears at 
0(5 4 ) - L-j is O(l) in 1/N C (see B). The rest of terms all remain subleading. 

The masses and the decay constants can be obtained from the two-point function of axial 
currents, which is most easily done by taking functional derivatives with respect to a^. We shall 
thus keep only terms quadratic in a M . The part of the effective action quadratic in the nonet fields 
finally reads 

S = \j d 4 x ([d^ a - b*)A ab (d^ b - 6j) - rB ab ^ b + rCabd^) , (7) 

where 6 M = 2/a M , cj) a (a = 0, 1, ...,8) are the fields with SU(3)r + l quantum numbers and 

A = I + AA , B = 2mB(D + AD) , C = AC . (8) 

" In the U(3)l®U (3)h chiral lagrangian, instead of constants the Li 's are functions of the 770 field. When we write 
Z/5 or Ls throughout the article we mean the value of those functions with the argument set to zero, £5(0), Lg(0). 
They are not to be confused with the constants in JJ, for which we reserve a superscript of [SU] L^ u \ 
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The mass matrix D is non-diagonal already at leading order: 

D n = D 22 = D 33 = 1 , D 4 4 = D 55 = D m = D 77 = l + ^, 

Dg8 = 1 +3 X ' D os = — , A)o = 1 + -^x - -— , (9) 

with m = (m u + m^)/2, x = (m s — m)/m. On the other hand, the matrix AA brings the next-to- 
leading order corrections to the kinetic term, 

AAn = AA 22 = AA 33 = 16^L 5 , 

AA 44 = AA 55 = AA m = AA 77 = 16-^-L 5 ^1 + - 

mB ( 2 
Avl 88 = 16-^L 5 \1 + -x 



AA 08 = j2-L 5 x , 

mB { 1 

A^ Q0 = 16— L 5 [l + -x\ , (10) 



and AD corrects the mass matrix, 



AD U = AD 22 = AD 33 = 32— L 8 , 



.mB 

T 2 



AD M = AD 55 = AD m = AD 77 = 32^L 8 ( I + .,■ - - r 

mB f 4 2 r 
AD 88 = 32-^ L 8 (l + -x + -x- 



AD 0S = j^L 8 x(2 + x) + V2v 3 ix , 

AD m = (l + ^x + ^x 2 ) - 2v 31 (3 + x) . (11) 



32^2 mB 
1 

mB , 2 I 

7 

AC comes entirely from the ?7a(1) anomaly: 

AQx) = !«50 • (12) 

From the correlator of two axial currents we can read off the physical masses and the decay 
constants, the former from the location of the poles, the latter from their residues. Following the 
steps of ref . @ , the correlator is the second derivative of (0) with respect to when evaluated for 

that minimises the effective action; it boils down to a simultaneous diagonalisation of A and B 
for the masses, 

A = I F , B = F Jf M 2 F , (13) 
whereas the decay constants read (P is a diagonal index, a is non-diagonal), 

f F „ = f((F^(A+l)) pa (14) 
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Consistently, we perform a perturbative diagonalization, order by order in the expansion. First, 
we diagonalize to leading order and only D has a non-diagonal part; it is already diagonal in n and 
K ( we shall neglect the isospin breaking ), 



2mB = Ml 

2mB (l + | ) = Mi , (15) 



only the and 8 components do mix, 

2mB FqDFI = M 2 , (16) 

with 

/cos^o -smflcA (17) 
V sin #o cos #o / 

From (15) we can fix mB and x unambiguously. Upon diagonalization of ( ]i~6| ) we fix V02 so as to 
describe correctly M 2 ,. This implies 

6»o ~ -20°, v 02 ^ -0.22 GeV 2 . (18) 

In this way the prediction for M 2 ~ (549 MeV) 2 which is only a 10% off. This is a known feature 
of the leading 1/N C approximation: it was proven in ref. that the ratio M 2 /M 2 , always comes 
out too small; the next-to-leading corrections are thus needed to reconcile the 1/N C expansion with 
the observed masses pX| ] . 

Let us emphasize that we had two masses, M 2 and M 2 ,, to fit and just one parameter, V02, to 
tune. In principle, we had the choice to fit either one to its experimental value. We could have, 
instead, fixed the value of vq2 in order to describe correctly M 2 , obtaining a prediction for M 2 , . 
This choice would give, though, a poor prediction for M 2 , since it is the rj that gets the bulk of 
the contribution from vq2 being, therefore, more sensitive to it. 

The next-order expressions for the masses can also be found in ||. The A A corrections are 
absorbed into a wave-function renormalization and the it and K are still diagonal 

Ml = 2mB (l + 16^(2L 8 - L 5 ] 

M| = 2m5 (l + |) (l + 16^(2L 8 -L 5 )(l + |)) , (19) 

whereas the and 8 components remain to be diagonalized and can be written as 
, AM 2 K -Ml 1 4 2 2 

m 2 8 = =^(M 2 K -Ml)(l + A M -A N ) , 
2 = l(2M 2 K + Ml)(l - 2A N ) + l(M 2 K - Ml) A M -3v 2 , (20) 



where 



m 00 j\ ■ — »/ ■ g 



8 , , ,0 , ,0, , . . _ 12 



Am = {M 2 K - Ml)(2L 8 - L 5 ) , A w = 3v 31 - ^v 02 L 5 . (21 ) 
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Note the explicit violation of the Kaplan-Manohar symmetry [11] in SU(3)l (8) SU(3)r - usually 



stated as the fact that in the expressions for the masses Lg and L5 always come through the 
combination (2L^ — L^) -m.U (3) chiral perturbation theory which is apparent in A^: this correction 
brings about a different dependence, on L5 only, due to the presence of V02, i-e., the Ua(1) anomaly. 

Let us now proceed to diagonalise to next-to-leading order. In order to preserve the form of the 
above partial result, we write the transformation matrix as 



cos 9 — sin 9 \ / 1 + 



2 2 

sin 9 cos 9 J V A408 1 + 



(22) 



where we have explicitly separated the rotation that diagonalizes the kinetic term. Note that 9 
also gets corrections with respect to 9q. For the physical masses we find, 



Ml = (1 ~ f )M| + \Ml - V9 + f + yl {Ml - Ml) + ( 1 - , ,„ ' I ' , ) 1 Ml - Mi ) Xu 




+ [ ~{2Ml + Ml)+ n lni \ —(3(2^ - 3M^) + y(-2M£ - ikffi ) ) A A . (2:!) 




where 



and 



3V9 + 27/ + 2/ 



V =H ^ (24) 



+ M*, = --(-3 + y)M 2 K + -yM^ + 2(M^ - M^)A M - -{2M Z K + M£)A N ; (25) 
the angle 9 being, 



tan 2, = M + ^ Am _ Ajv _ 1 ^±^A„) . (26, 
l + y \ l + y l + y Mfc - M£ J 

Let us now turn to the numerical exploitation of these results. The global counting of parameters 
versus data goes as follows: we are left with y (which is proportional to ^02), Am and A at to fix 
M^ and M^,. Thus, no prediction can be made. Yet, if we fix a rotation angle, then there is a fit 
for all three parameters. The results for angles which are close to the zeroth order #0 are displayed 
in the following table: 



9 


A M 


A N 


VQ2 


-18° 


0.113 


0.270 


30.8 


-20° 


0.156 


0.220 


29.3 


-22° 


0.203 


0.178 


27.9 


-24° 


0.254 


0.143 


26.5 



The best fit, if U(3) chiral perturbation theory is to make sense, is expected to be near the 
zeroth order result. One could play with different convergence criteria to argue in favor of a given 
angle. A nice feature of our table is that Am and A^v move in different directions as we change the 
mixing angle. This leads to a minimization of corrections around 9 ~ —20°. Moreover, one may 
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select the optimal 9 where the masses of r\ and rf receive the smallest correction in the sense that 
corrections proportional to Am and An tend to cancel. Interestingly enough, we find it to be very 
close to —20° as well. 

Of course none of these arguments is a substitute for a global fit to data. However, our results 
agree in hinting at a mixing angle close to —20° > 9 > —21°, which correspond to values of 
and A^v compatible with those given in § (A M ~ 0.18, Aiv ~ 0.24). We stress that since from 
this point of view the angle 9 is the only variable, the values of Am and A^v are correlated. 

At this point let us comment on the Gell-Mann-Okubo formula, with corrections from chiral 
perturbation theory. The violation of the original formula is related to Am, 

AM K -Ml- 3M 2 - 3 sin 2 9 (m 2 , - M 2 ) = -4(M| - M 2 )A M . (27) 



The leading order equation, thus, corresponds to setting the r.h.s. to zero. At this order, 9 ~ — 10 



o 



in order to fulfil the equation. This is not the procedure we chose in Eq.18, as we took M 2 , to fix 
9 ~ —20° at leading order. Eq. (|27]) thus provides a different way of fitting the mixing, which is 
strongly corrected by Am- Of course, at next-order, the best fit at e.g. Am ~ 0.156 is 9 ~ -20° 
as well. 

We may also use our analytical results to get the form of the decay constants for the whole 
nonet. The diagonal elements of relevance are given by the combination 

/cosfl -sin0\/l + ^4 M ^f* \fcos9 sm0\ . , 

Iphysical -J[ s - m0 cos Q ) [ AAca j_ + + ) { _ gin q cos q) ' ^> 

This expression yields 

7 
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U = f ( 1 + 4 tI m ^ 



Sk = f ( 1 + 



W*,2 M K~ M l 9 + ^ \ , 3 A l + y 



fv = f 1 + 4^ \M' K - " /nin \ 5 + o«50 1 + 



P \ 3 ^9 + 2y + y* J 8 ^ + 2 V + V 2 ) ) ' 

^/( 1 + ^(MUM_^) + 3„ 5o(l __^_)) , 29) 

Note that in the exact SU(3)r+l limit of quarks degenerate in mass there is no mixing and one 
finds 



/„ - / 1 + 4 



L 5 4MI-M 2 ' 
P 3 



f t h i ^ 2M K + Ml 3 \ 

frf -> f I 1 + 4 ^2 § + 4 : V50 J ' ( 30 ) 



which is the y — > — oo limit of (|29|), as expected. 

Let us compare these results with data. We can take f n and fx to fix / and L5, (f n = 92 A 
MeV, j K = 1.223/,,) 

/ = 90.8MeV , L 5 = 2.0 10~ 3 . (31) 
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This is the same value as obtained for L<\? U ^ r in [|j for it is extracted from the same source. In 
our approximation the issue of the running of the L5, L$ cannot be addressed because the loop 
corrections are dropped altogether. As a matter of principle one expects L5 ~ L^ U ^ because 
integrating out the rf does not give any contribution to l\? U ^ , L^ U ^ at tree level, except for Ly St/ ' . 

We are left with f n and /„> to be described by y (again, proportional to V02) and U50. Assuming 
that convergence criteria fix a rotation angle around 6 ~ —20°, this leaves room for a fit of V50 and 
a prediction. The values of 77 and rj decay constants are, though, poorly determined. Indeed, the 
experimental values of f n and f v i || have large errors: 

0.943 < ^ < 1.091 , (32) 
0.912 < 4^ < 1-015 . 

Let us present the output of our formulae as a function of the mixing angle 6. At each given 9, 
the experimental values of fy- allow for a range of values of the parameter v 50 and of -2- . We again 
display the results for angles which are close to the zeroth order 0q. 



e 


— ^50 min/max 


Y~ min/max 


-20° 


2.293/0.601 


-0.218/0.902 


-22° 


1.775/0.365 


0.171/1.080 


-24° 


1.394/0.197 


0.457/1.205 



A word of caution should be added as regards the use of the /„/ value, as made in this paper. 



As pointed out by Shore and Veneziano [12|, from the decay rate of rf — ► 77 one cannot obtain 



(5) 

the value of because the singlet axial current is not conserved in the chiral limit, due to 

(5) 

the £7a(1) anomaly, which makes from (0\Afj, \n (k)) = ik^frf depend on a subtraction point 
and thus not be observable. Of course this dependence is 1/N C suppressed because in the large- iV c 
limit the anomaly is a subleading effect. The relation between our /„/ and what is measured in 



v' ~^ 77 an d given in ||, [13] can be worked out and ammounts to include one further term in 
the effective action coupled to an external electromagnetic source 0. This would still bring in 
a new unkown constant and we would face the problem of having one more constant to fit than 
measured constants available. Nevertheless, we overcome this shortcoming by using the criterium 
of minimum sensitivity that minimizes the size of the corrections to a given order and conclude that 
the combined expansion @, within the framework of Ul(3) <8> Ur(3) chiral perturbation theory, 
is able to accommodate very naturally the observed values of masses and decay constants /„-, fx- 
The big uncertainty in makes the method inconclusive for (which, strictly speaking, is not 
experimentally known), as can be seen from the last table. If data from 7/ — > 77 were more precise, 
we could with our method pursue to determine the new contant that is involved in the process, 
and how well it would accommodate in the framework of the combined expansion. 

We finish by quoting the values of all the parameters in one batch, as functions of the mixing 
angle in the range 20° < 6 < 24°: 

0.15 < A M < 0.26, (33) 



8 



0.980 < 
18.3 < 
-4.7 < 
26 < 
0.14 < 
1.35 10 -3 < 
-0.164 < 



We should like to add that after this work was finished there appeared the document of ref. 
|i"4f in which the author reports some results, as yet unpublished, on the same issue that concerns 
our study. 
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